Chandra observations of the massive star-forming region S106 
X-ray emission from tlie embedded massive protostellar object IRS 4 
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Abstract. We present Chandra observations of the massive star-forming region S106, a prominent Hn region in Cygnus, 
associated with an extended molecular cloud and a young cluster. The nebula is excited by a single young massive star located 
at the center of the molecular cloud and the embedded cluster. The prominence of the cluster in the Chandra observation 
presented here confirms its youth and allows some of its members to be studied in more detail. We detect X-ray emission from 
the young massive central source S106 IRS 4, the deeply embedded central object which drives the bipolar nebula with a mass 
loss rate approximately 1-2 orders of magnitude higher than main sequence stars of comparable luminosity. Still, on the basis 
of its wind momentum flux the X-ray luminosity of SI 06 IRS 4 is comparable to the values observed in more evolved (main 
sequence and giant) massive stars, suggesting that the same process which is responsible for the observed X-ray emission from 
older massive stars is already at work at these early stages. 

Key words. Stars: activity - Stars: formation - Stars: pre-main sequence - Stars: individual: SI06 IRS 4 - X-rays: stars 



1. Introduction 

The study of regions of recent star formation provides the ob- 
servational data necessary to address a number of important 
astrophysical questions, such as the physics of the star forma- 
tion process, what determines its efficiency, the characteriza- 
tion of the initial mass function, the stellar chemical evolution 
and the galactic chemical recycling process; in particular, mas- 
sive stars, with their strong winds and short lifetime have a 
strong impact on the circumstellar environment. The high den- 
sity of absorbing material makes optical observations poorly 
suited for the observation of the earliest stages of star forma- 
tion, while the infrared (IR) bands J, H and K have proven to 
be an excellent tool for unveiling embedded star formation re- 
gions due to the reduced impact of dust extinction at longer 
wavelengths. 

A role similar to IR observations is currently being played 
by X-ray observations, to which dust is as transparent. Given 
that X-ray luminosity of stars evolves strongly with age, young 
stars are very active X-ray sources, so that they can be easily 
selected from field stars. X-ray surveys of young associations 
and clusters also provide valuable information on the magnetic 
activity of young stars and on their accretion rates, as well as ef- 
fectively supplement other approaches to membership determi- 
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nation in these regions (see e.g. lpeigelson & Montmerlel 19991 
although the field is evolving rapidly). In comparison with the 
low-mass domain X-ray observations of massive star-forming 
regions and massive young stellar objects (YSOs) have been 
been rare, with very few detecti ons of X-ray em ission from 
massive YSOs reported to date jKohno et alJ2002h . 

Recentlv'Le Du igou & Kn odlseder' ( l2002l) (LDK02, here- 
after) have published a study (based on the 2MASS Point 
Source Catalogue - PSC) of 22 young open clusters in the 
Cygnus region, of which 12 objects are recently discovered 
clusters and 3 are new clusters candidates. All of them ap- 
pear to contain a significant number of massive stars. We have 
searched the Chandra and XMM-Newton archives for observa- 
tions of these regions, and found a Chandra medium sensitiv- 
ity observation (45 ks) in the direction of one of the recently 
discovered clusters, which corresponds to the well known re- 
gion of star formation S106 (referred to as "C102" in LDK02). 
Note that there are no Chandra or XMM-Newton observations 
to date in the direction of any of the other 21 young clusters. 

The Hn region SI 06 is a prominent bipolar emission neb- 
ula associated with an extended molecular cloud. The nebula 
appears to b e excited by a sin gle massive star at its center, 
SI 06 IRS 4 dGehrz et alJll982 '). also known as SI 06 IR. The 
star is surrounded by a a compact region of radio emission indi- 
cating that this is a newly formed star still embedded within it s 
natal dust cocoon JSimon et alJl983tlCrowther & Contil2003l) . 



IStaude e t al.l (11^)82 ) have used U Vail photometry and broad- 
band polarimetry of field stars to determine the foreground ex- 
tinction and polarization, resulting in a distance estimate for 
the molecular complex of 600 + 100 pc {DM = 9). 

Th e SI 06 cluster was first studied bv 'Hodap p & Ravnen 
11199 11) in the IR, who counted ~ 160 cluster members of which 
half are within a radius of 1 .0 arcmin from the cluster centre - 
which they position at about 30 arcsec west and 15 arcsec north 
of S 106 IRS 4. They derive a cluster mass of ~ 140 Mq and es- 
timate that star formation in this cluster has been ongoing for 
the past 1-2 million years. 

In their study, LDK02 position the cluster center at 
a, 5=20:27:25, 37:22:48 (all coordinates are in the J2000 
system), which is about 21 arcsec west of the position of 
S106 IRS 4 in the 2MASS PSC (the coordinate of S 106 IRS 4 
being a, 5=20:27:26.76, 37:22:47.9). From cumulative star 
counts they derive a 90% population radius TJgo = 1-1 arcmin, 
even though they caution that this value may only correspond 
to the core radius of the structure. They estimate the cluster to 
contain 20 + 5 OB star and have a total mass in the range of 
400-600Mo. Schneider et al. (2002) identify the cluster as the 
primary star formation site in S106, coincident with the promi- 
nent peak in CO emission. They also detected a second weaker 
CO peak located 5 arcmin south of SI 06 IRS 4, which harbors 
a small IR cluster and nebulosity. They interpret the nebula as 
a signpost for a secondary site of star formation in SI 06 and 
refer to it as "S106 south". 

In this article we present Chandra observations of the SI 06 
and SI 06 south clusters and report detection of SI 06 IRS 4 in 
the X-ray. The paper is organized as follow: the observations 
and data analysis are presented in Sect.|2| The results are sum- 
marised in Sect. Inland discussed in Sec.|4l 



2. Observations 

The X-ray observations discussed in this paper were obtained 
with the Chandra observatory. The 45 ks ACIS-I observation of 
S106 was taken starting on November 3 2001 at 00:08:14 UT. 
The Principal Investigator for these observations is Y. Maeda 
and the observation target is the protostellar object S 106 IRS 4. 
The data were retrieved from the public data archive, with no 
re-processing done on the archival data. 

We performed the source detection on the event list, us- 
ing the Wavelet Transform detection algorithm developed 
at Palermo Astronomical Observatory (Pwdetect, available 
at http://oapa.astropa.unipa.it/progetti_ricerca/PWDetecti. The 
threshold for source detection was taken as to ensure a maxi- 
mum of one spurious source per field. 

For the brightest sources for which a spectral and timing 
analysis was carried out source and background regions were 
defined in ds9, and light curves and spectra were extracted 
from the photon list using CIAO V. 2.2.1 threads, which were 
also used for the generation of the relative response matrices. 
Spectral analysis was performed in xspec. 

Fig.^shows the field centered on S 106 as seen in the ACIS- 
I camera (for various energy ranges) and in the /T-band of the 



/MASS survey . Only the area withm 1.1 arcmm from the cen- 
tre is displayed for clarity, although we applied our analysis 
to an area within 4.4 arcmin from the cluster centre. A radius 
of 4.4 arcmin was chosen with the aim of studying the clus- 
ter with a minimal contamination from foreground and back- 
ground sources in the absence of a membership list. The chosen 
radius corresponds to 4 times the SI 06 core radius as estimated 
by LDK02 and therefore should contain the large majority of 
the cluster members. At the same time, it does not extends to 
regions too far off-axis, where the sensitivity to X-ray sources 
would be lower (because of vignetting and because the Chandra 
PSF widens). 

A clustering of X-ray sources is visible in the ACIS-I image 
obtained for the entire energy range (AE = 0.2 - 10.0 keV), al- 
though the cluster is 'cut' by the ACIS CCD gaps, which deter- 
mine a reduced sensitivity region. The 2MASS ^-band image 
is dominated by the young massive stellar object SI 06 IRS4 
(a bright IR source) near the center and by the double-lobed 
infrared nebula around it. SI 02 IRS4 is thought to be the en- 
gine which po wers the doubl e-lobed nebula and the associated 
H n region (Ge hrz et al.ll982l) . The X-ray sources detected with 
Pwdetect are overlaid on the 2MASS image. 

As discussed below, most X-ray sources (and in partic- 
ular SI 06 IRS 4) have an IR counterpart. A number of X- 
ray sources close to SI 06 IRS 4 have no counterpart unre- 
solved from SI 06 IRS 4 in the 2MASS image, while north of 
SI 06 IRS 4 an example of a double X-ray source with an un- 
resolved IR counterpart can be seen. The cluster center coordi- 
nates as given by LDK02 (the center of the big circle in Fig.[0 
lies 0.3 arcmin off-axis in the Chandra observation. 

In Fig.[2the ACIS-I images of the field centered on S106 
for a soft {JS.E - 0.2 -1.0 ke V) and a hard energy band (AE — 
1.0 - 7.5 keV) are also shown. The great majority of the X- 
ray sources in the direction of SI 06 are not visible for energies 
lower the 1 .0 keV, while are clearly visible in the hard band, as 
expected on the basis of the high extinction toward the cluster 
(see below). 

In Table^we list the coordinates for all the X-ray sources 
detected with Pwdetect within a radius of 4.4 arcmin from the 
center of S106. A total of 87 X-ray sources are detected, of 
which 45 are within the SI 06 core radius. The X-ray source 
density within the 4.4 arcmin radius is 1.4 arcmin"^, 10 times 
higher than in the rest of the ACIS-I field^. The X-ray source 
density within the 1.1 arcmin core radius is 11.8 arcmin"^. 
Most of the detected sources are rather weak, with a count 
rate below 1 ct ks"'. The faintest source is source 9 with 
a count rate of 0.16 ks ', which corresponds^ to a flux of 
1.6 X 10 '^ erg cm"^ s"'and an intrinsic (unabsorbed) flux of 
4.2 X 10"'^ erg cm"^ s"'. Only 7 sources (sources 22, 30, 32 
39, 60, 68 and 72) have a count rate just above 2 ct ks"'. 



' obtained from the 2MASS Quicklook Image Services of the 
NASA/IPAC Infrared Science Archive. 

^ excluding the sources and the area of "S106 south" ~ see next 
section. 

' Assuming a plasma temperature of ^r = 2.16 keV, a metallicity 
of Z = 0.2 Zq and N{ii) as derived from IR data for this source. See 
Sect. 3.2. 




Fig. 1. The field centered on S 106 in the ACIS-I camera for the full energy range of the detectors (top left), in the2MASS /T-band 
(top right), in the ACIS-I camera for the energy range A/s = 0.2 - 1.0 keV (bottom left) and in the ACIS-I camera for the energy 
range AE = 1.0 - 7.5 keV (bottom right). The X-ray sources detected in the Chandra image have been overlaid on the 2MASS 
image and are represented by small circles. The core size of SI 06 is shown as a large circle (radius 1.1. arcmin). In the top left 
image the lines indicate the position of the CCD gaps in the ACIS camera. 



For the sources with an IR counterpart the values of their 
magnitude in the J, H and K bands are also given in Table ^ 
together with the radial distance from the possible IR counter- 
part. An X-ray source was considered to have a possible IR 
counterpart if an IR source is present in the 2MASS PSC cat- 
alogue within a radius of 3 arcsec from the X-ray position. Of 
the 87 sources within the 4.4 arcmin radius 64 have a possible 
IR counterpart. Source 82 has two possible IR counterparts, 
while sources 55 and 58 have the same IR counterpart. All the 
7 brightest X-ray sources have an IR counterpart. 

The X-ray counterpart of S 1 06 IRS 4 is source 50 in our list. 
It is a weak source with a count rate of 0.3 + 0. 1 1 ct ks"' ; its X- 
ray position is less than 0.5 arcsec away from the coordinates of 
S106 IRS 4 in the 2MASS PSC, which agrees we l l with the co- 
ordinate for S 106 IRS 4 given in lSchneider et alJ (|2002') (a, 5^ 
20:21:26.1 A, 37:22:48), but are more than 12 arcsec away from 
the position given for this object in the Simbad database, w hich 
corresponds to the coordinates given in lOehrz et alJ 1119821) . 



To quantify the visual impression that the great majority of 
the X-ray sources detected in the direction of SI 06 have hard 
spectra, (and thus are not visible below 1 .0 keV) we have run 
PwDETECT also on the energy-filtered ACIS-I event lists. Only 
8 sources (flagged with a 's' in Table ^ are detected by the 
algorithm in the soft energy band {AE - 0.2 -1.0 keV) within 
a radius of 4.4 arcmin from the center of S106. All of these 
8 sources have an IR counterpart, and only one is within the 
S106 core radius. In the total ACIS-I field 27 out 125 sources 
are detected in the soft band, implying that, within 4.4 arcmin 
from the center of S106, the fraction of soft sources (9%) is 
more than 6 times lower then in the rest of the ACIS-I field 
(60%)! 

LDK02 derive, for the stars belonging to S106, extinction 
coefficients in the K band ranging from A/f = 0.5 to Ak - 4.0. 
Given the conversion factor cRieke & Lebofskv-1985.) : 



O.llAv 



(1) 



^ excluding the sources in "S106 south" - see next section. 



and the relation (see e.g. lCoxIzUUUD : 
NiU)/Av = 1.9 X 10^' cm-2 mag-' 



(2) 



the absorbing column density for members of S106 must be 
N(}i) S 8.5 X 10^' cm"^. The fact that the majority of the 
sources are not visible in the ACIS camera for energies be- 
low 1.0 keV is consistent with such high absorbing column 
densities. Indeed, if we assume as a typical spectrum for the 
sources listed in Tableman absorbed Raymond-Smith plasma 
model with A^(H) = 2.0 x 10^' cm"^, plasma temperature 
kT = 1.5 keV and metal abundance Z - 0.2 Zq then, using 
the piMMS simulator at Heasarc, the expected ratio of the count 
rate between the soft and the total band is 0.15. This estimate 
is not very sensitive to the assumed X-ray temperature in the 
range 1- 3 keV. 

2.1. S1 06 south 

As mentioned in the Introduction. lSchneider et alJ 1120021) iden- 
tify the cluster SI 06 as a primary star formation site coincident 
with the prominent peak in CO emission. They also detect a 
second weaker CO peak located 5 arcmin south of SI 06 IRS 4 
which harbors a small IR cluster and nebulosity and refer to it 
as "SI 06 south". 

Fig. 121 shows the field centered around SI 06 south in the 
ACIS-camera and in the 2MASS K band. An increase in source 
density in the direction of S 106 south is visible in the Chandra 
image. The approximate center of this group of X-ray sources 
is at a, 6 = 20:27:23.2, 37:17:33 (5 arcmin off-axis in the 
Chandra data). In Table |2lthe coordinates and 2MASS IR pho- 
tometry for all X-ray sources within a radius of 0.8 arcmin from 
this position are listed. A total of 6 sources have been detected, 
making the source density within this area 20 times higher than 
in the rest of the field (excluding the area of 4.4 arcmin radius 
from the center of S 106). 

All of the 6 sources have an IR counterpart within a radial 
distance of 3 arcsec, and only one of them is visible in the soft 
band. 

3. Results 

3.1. Color-magnitude and color-color diagrams for the 
IR counterparts to X-ray sources 

To compare the population of the X-ray sources detected in the 
Chandra exposure with the population of SI 06 as defined by 
LDK02 we have constructed a K \s. J - K color magnitude 
diagram (Fig.O for the X-ray sources listed in Tables^andlJl 
which have an IR counterpart in the 2MASS PSC. 

In the diagram 2 Myr isochrones are shown for 4 values of 
extinction (Ak = 0, 1, 2, 3). Reddening vectors have been com- 
puted following the same relation adopted by LDK02, Ak = 
Rk X E{J - K) where E(J - K) is the colour excess (E(J - K) - 
(J-K)-(J-K)o) an d Rk = 0.66 iJRieke & Lebofskvll985t) . The 
isochrones are from lSiess et al.l ( 12000 1). for a metal abundance 
Z - 0.02 plus overshooting, shifted to the estimated distance of 
S106 (600 pc). Various simplifying assumptions apply to these 
models, e.g. they include neither rotation nor accretion, and we 



caution that evolutionary models tor pre-main sequence stars 
are not yet well established. Stassun et al. (2004 ) recently re- 
ported the discovery of a double-lined, spectroscopic, eclipsing 
binary in the Orion star-forming region, with measured masses 
of 1.01 Mq for the primary and 0.73 M0 for the secondary. 
They used their measurements of the fundamental stellar prop- 
erties for both components to test the predictions of pre-main 
sequence stellar evolutionary t racks. None of th e models they 
examined (including the one bv lSiess et alJ200f)l) correctly pre- 
dicts the masses of the two components simultaneously. 

In Fig. |3] seven sources have been highlighted using a tri- 
angle (sources 16, 27, 35, 67, 68, 80 and 84). These sources 
fall on the blue side of the 2 Myr isochrone for zero extinction 
and therefore are inconsistent with the estimated age and dis- 
tance of S106. There is evidence that these maybe foreground 
sources. They do not meet the membership criteria used by 
LDK02 (their extinction coefficient being too low) and all of 
them but source 16 are detected in the Chandra soft band, con- 
sistent with their being less absorbed than the others. Note, 
however, that sources 30 and 87, that have been flagged as 
"soft", are not inconsistent with a 2 Myr isochrone and the 
value of their extinction coefficient as computed by LDK02 is 
consistent with their membership criteria. 

The sources in our sample appear very scattered in the 
color-magnitude diagram, implying highly variable absorption. 
LDK02 find the same large scatter for their IR sample of 
sources belonging to SI 06 (according to their selection crite- 
ria), and suggest that, while this scatter is partly due to large 
and variable absorption, it is probably also an effect of the un- 
reliability of the 2MASS photometric data in this very crowded 
area''. From the color-magnitude diagram it appears that the 
X-ray selected sample contains a number of heavily absorbed 
massive stars. The brightest of all, with K - 5.9, is the mas- 
sive young stellar object S106 IRS 4 at the center of the S106 
complex. 

The color-magnitude diagram also provides an indication 
of the completeness limit of our X-ray survey, which is at 
around K ^ 12, as we have verified by plotting differential 
star counts versus apparent magnitude. For a 2 Myr isochrone 
K = 12 corresponds to stars of mass M ^ 0.5 Mq for low ex- 
tinction [Ak < 0.1) and M ^ 1.2 Mq for extinction coefficients 
Ak > 1.0 

The J - H versus H - K color-color diagram in Fig. [fal- 
lows normally reddened stars to be discriminated from star with 
IR excess (indicative of warm circumstellar dust in addition to 
a reddened photosphere). The intrinsic colours for stars on a 
2 Myr isochrone are indicated by the solid line, whereas the 
dotted line yields the locus of dereddened colo urs of classi- 
cal T Tauri stars according to lMever et al.l (Il997l) . The dashed 
lines define the region of normal reddening using the redden- 
ing law from Rieke & Lebofskv ( 1985) (the normal reddening 
region of main sequence stars would be very similar to the one 
of the stars on the 2 Myr isochrone). The scatter of ffie X-ray 
sources outside this region is significant and it is probably an 
effect of the unreliability of the 2MASS photometric data in 
this crowded area. Nevertheless in the diagram we have identi- 



The 2MASS images resolution is 1 arc-sec/pixel. 
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Fig. 2. The field centered on S106 south in the ACIS-I camera (left) and in the 2MASS /T-band image (right). The X-ray sources 
detected in the Chandra image have been overlaid on the 2MASS image and are represented by small circles with size propor- 
tional to the X-ray count rate. The selection area that we have used for S106 south is shown as a large circle (radius 0.8 arcmin). 
In the ACIS image the lines indicate the position of the CCD gaps. 



fied 10 sources which appear to have significant IR excess (i.e. 
lie to the right of the normal reddening region) and therefore 
could be embedded young objects. The presence of an IR ex- 
cess for source 50 (S106 IRS 4) is akeady known from the lit- 
eratur e (IFelli et al..l984.: .Gehrz et al,..1982; Hodapp & Rayner 
Il99lh . 

In the color-magnitude diagram sources 29, 50 
(S106 IRS 4) and 62 appear to be massive objects (with 
masses above 7 Mq). This estimate does not take into account 
the presence of an IR excess which may lead to overestimate 
the star mass. This is not the case however for S106 IRS 4, for 
which estimates in the literature independently indicate a mass 
greater than 15 Mq JFelliet alJl984 . 



3.2. X-ray luminosity as a function of stellar mass 

To convert detector count rates to intrinsic X-ray luminos- 
ity we have derived, using the pimms software at Heasarc, a 
conversion factor between (unabsorbed) source flux and ob- 
served count rate, taking into account the absorption. The ab- 
sorbing column density for each star was derived by first es- 
timating Ak and then deriving A^(H) using Eqs. Q and (|2}. 
The value of Ak was determined assuming that all sources in 
S 106 and S 106 south would intrinsically lie on the theoretical 2 
Myr isochrone and then measuring the displacement along the 
reddening vector No extinction can obviously be derived for 
sources 16, 27, 35, 67, 68, 80 and 84, which lie on the blue side 
of the isochrone. 



Because of the shape of the 2 Myr isochrone a degeneracy 
is present for mass ranges 2.2-7.0Mq and 0.2-0.3 Mq; lack- 
ing additional information (e.g. spectral types or optical magni- 
tudes) this degeneracy cannot be resolved. We thus, somewhat 
arbitrarily, decided to associate the mass derived through the 
rightmost intersection with the isochrone for the mass range 
2.2-7.0Mq, and redefined the isochrone for the mass range 
0.2-0.3 Mq using a tangent along the reddening vector at M = 
0.2 Mq (shown in Fig.|3]as a thin line). 

As a typical X-ray spectrum for our sample we assumed 
a one-temperature plasma model with kT - 2.16 keV* and 
a coronal metal abundance Z = 0.2Zo, deriving, for A^(H) - 
1 .0x10^^ cm"^ and a count rate of 1 ks"' , a conversion factor of 
^x intr - 2.7 X lO^'^^erg cm"^ s"'. The conversion factor scales 
almost hnearly with N(U) (for N(ii) = [0.5 - 6.0] x lO^^ cm^^, 
the range in which all our objects fall, see below), so that 



Fx inti- = 2.7 X 10"'^ X A^(H) x ct [erg cm"^ s"'] 



(3) 



where ct are the source counts in units of ks"' and A^(H) is 
the absorbing column density in units of 10~^ cm -. The error 
introduced in the flux estimate by the linear approximation is 
below 30%. 



^ The value is the same as the one assumed by iFlaccomio et alJ 
1(2003) for their analysis of stars in the ONC, to which we will be 
comparing our results; this is similar to the average plasma tempera- 
ture of kT = 2.4 keV derived for the brightest sources in Table 3 (next 
section), excluding source 30 (because of the large uncertainty) and 
source 68 (because likely a foreground source) 
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Fig. 3. Color magnitude diagram for the IR counterparts to 
the X-ray sources detected within 4.4 arcmin from the center 
of S106 (black circles) and within S106 south (gray squares). 
Theoretical 2 Myr isochrones are also shown for different ex- 
tinction coefficients (from left to right Ak - 0, 1,2,3). The 
dashed lines indicate the reddening vector for stars of constant 
mass. Points with a triangle indicate sources inconsistent with 
the theoretical 2 Myr isochrone (see text for details). 



Fig. 4. Color-color diagram for the IR counterparts to the X-ray 
sources detected within 4.4 arcmin from the center of S 1 06 (cir- 
cular black points) and within S106 south (square gray points). 
The solid line at lower-left is the theoretical 2 Myr isochrone. 
The dotted line yields the locus of dereddened colours of clas- 
sical T Tauri stars according to Meyer et al. (1997). The dashed 
lines mark the reddening bands. The most heavily reddened 
stars are identified. 



In all but three cases we derived a value of A^ correspond- 
ing to an absorbing column density within the range [0.5 -6.0] x 
10^^ cm"^. For sources 38, 82-a and 87 we derived values of 
Ak < 0.07 (corresponding to A^(H) < 0.1 x 10^^ cm"^), signifi- 
cantly below the absorbing column densities derived for all the 
other sources in our sample (which has a median of A^ = 1.1); 
we have therefore excluded these sources from the rest of the 
analysis. Having determined for a star its intrinsic magnitude in 
K, as described above, an estimate of its mass is readily avail- 
able by using the 2 Myr isochrone. 

The derived X-ray luminosity as a function of the stellar 
mass is shown in Fig.|5j 61 sources are included in the final 
samples, i.e. all sources in Tables^andQwith the exception of 
sources 16, 27, 35, 38, 67, 68, 80, 82-a, 84 and 87. The median 
X-ray luminosity for stars in different mass intervals is also 
shown , using the same mass bins adopted by Iplaccomio et alJ 
1I2OO3I) in their analysis of the X-ray activity indicators for 
members of the Orion Nebula Cluster (ONC). The median X- 
ray luminosities are summarised in Table |3] We note that for 
masses lower than 0.5 Mq the derived median X-ray luminosi- 
ties are systematically higher than the corresponding ONC stars 
while the reverse is true for M > 0.5 Mq. The values at lower 
masses (M < 0.5 Mq), which corresponds to lower fluxes, are 
likely biased to higher X-ray luminosities due to lack of a par- 
ent sample of cluster m embers. In their analysis of the ONC 
Iplaccomio et alJ 1120031) made use of an optical sample of clus- 



ter members in order to place upper limits on X-ray luminosi- 
ties. Lacking such a sample for the S106 clusters we are not 
able to place upper limits. The fact that for M > 0.5 Mq the 
values of median Lx for our sample are systematically lower 
than for the ONC could be an indication of some physical dif- 
ference between the SI 06 region and the ONC, such as a dif- 
ferent fraction of binaries. Nevertheless the statistic is too small 
and the uncertainties too large to attach any significance to this 
difference. 

A trend of increasing X-ray luminosities with star mass is 
clearly visible in Fig.|5] Two areas, one for 0.17-0.46 Mq and 
one for 2.9-9.4 Mq, are devoid of points. This is, partly, an 
effect of the way we have chosen to circumvent the degener- 
acy in value of A^ for some mass ranges (see above), partly, it 
is a feature of the data themselves. We note also that lacking 
optical/IR spectral studies of the members of the SI 06 clusters 
we do not know neither the spectral energy distribution nor the 
spectral characteristics of the stars in our sample and are there- 
fore unable to classify them. 

In Fig.lsjthe star with a mass estimate of 102 Mq is source 
50, corresponding to the central exciting star SI 06 IRS 4. The 
value of 102 Mq is very probably overestimated, due to the 
fact that the method we have used to derive a mass value for 
each star does not take into account the possible presence of 
an IR excess, wh ich in this case is known to be present (e.g. 
iFelli et al.lll984l and Fig. 0}. Indeed, the masses of all stars 
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Fig. 5. Scatter plot of intrinsic X-ray luminosities as a func- 
tion of mass for the sample of 61 sources defined in the text. 
The gray horizontal lines indicate the median Lx for the cor- 
responding mass interval. The abscissa of the empty-squares is 
the median mass for each subsample. The same mass intervals 
(in units of Mq) as in Flaccomio et a l. (2003) have been con- 
sidered: (0,0.16], (0.16,0.25], (0.25,0.5], (0.5, 1.0], (1.0,2.0], 
(2.0, 3.0], (3.0, 10] and (10, oo) 



in our sample for which we identified a significant IR excess 
maybe affected by such an error For SI 06 IRS 4 th e estimates 
in the literature indicate a mass greater than 15 Mq ( Fell i et alJ 
11984) and therefore the presence of an IR excess does not affect 
the placement of SI 06 IRS 4 in the highest mass bin. 

Table 3. Median value of the X-ray luminosity (second col- 
umn) for our sample of 61 sources in SI 06 and SI 06 south for 
different mass ranges. These values can be compared with the 
values fo r members of the ONC (fourth column), as derived by 
iFlaccomio et al. (2003) . The number of sources of our sample 
within each mass range is also given. 



Mass 

Mq 



log(Lx) 
log(ergs"') 



N. src. 



log(Lx) ONC 
log(ergs"') 



<0.16 


29.57 


0.16-0.25 


29.60 


0.25-0.50 


29.85 


0.50-1.00 


29.90 


1.00-2.00 


30.17 


2.00-3.00 


30.30 


3.00-10.0 


30.42 


> 10.0 


30.51 



2 


29.00 


3 


29.15 


3 


29.60 





30.05 


8 


30.45 





31.20 


1 


30.50 


4 


30.90 



J.J. ispectrai ana timing analysis or me ongntest A-ray 
sources 

A spectral and timing analysis was carried out for the seven 
brightest X-ray sources in the list of Tablen(sources 22, 30, 32, 
39, 60, 68, 72). Their light curves are shown in Fig.|6l Sources 
22 and 68 present large-amplitude variability, and source 22 is 
not visible (no photons are detected) for the first 10 ks of the 
Chandra exposure; the rise in the source luminosity at 10 ks 
is impulsive and it is followed by a decay lasting about 10 ks, 
with a behavior typical of stellar flares. The mass of source 22 
is estimated at M ^ 0.6 Mq. The time variability of source 68 
is remarkable, with the source counts suddenly increasing to 
more than 10 times the previous level in about 2 ks and then 
decreasing abruptly again to a level about 3 times higher than 
before the jump. We do not provide an estimate for the mass as 
its position in the color-magnitude diagram is inconsistent with 
a 2 Myr isochrone at 600 pc. Both in the X-ray image and in 
the 2MASS data source 68 has an apparent companion, source 
67, at about 4.7 arcsec. Both sources are visible in the X-ray 
soft band and their Ak derived with the approach described by 
LDK02 are very similar, consistent with the two sources being 
physically associated. 

The Kolmogorov-Smirnov (K-S) test, which measures the 
maximum deviation of the integral photon arrival times from a 
constant source model, was applied to the light curves, with the 
results summarised in Table|4] As expected sources 22 and 68 
have a negligible probability of being constant. The probability 
of constancy is also low for sources 30 and 39. 

The background subtracted spectra for the seven brightest 
sources are plotted in Fig.|6l together with the best-fit models. 
The ACIS-I spectra, rebinned to a minimum of 5 counts per 
energy bin, were fit with an absorbed one-temperature plasma 
model with me tal abundance 0.2 o f solar, a value typical for T 
Tauri stars (e.g. lFavata et alJ20()3l) . The best fit spectral param- 
eters are summarised in Table |3 The average coronal temper- 
ature for all sources but source 68 (likely a foreground source) 
is kT - 3.1 keV, excluding also source 30 because of the large 
uncertainty the average is kT - 2.4 keV. From the spectral fits 
intrinsic X-ray luminosities have been derived (also shown in 
Table|4}, assuming a source distance d = 600 pc. 

Excluding source 68, the values of absorbing column densi- 
ties derived for these sources are all higher than 5.0 x 10^' cm"^ 
(corresponding to Ak > 0.3), consistent with the sources being 
in a region of high extinction. Sources 30 and 32 have values 
of A'(H) around the lower limit of the typical absorbing column 
density for our sample of 61 sources, so that they could also 
be active foreground stars. Nevertheless their Ak is consistent 
with cluster membership, as is their X-ray luminosity. 

The best-fit values of A^(H) and Lx reported in Table |3 are 
consistent, within a factor of three for the X-ray luminosity and 
within 2cr for A^(H), with the ones computed using the pho- 
tometrically derived A^(H) and assuming a typical spectrum. 
This provides a consistency check on the method used to de- 
rive A^(H) and Lx for all other sources in our sample of 61 for 
which a spectral analysis was not possible. 
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Fig. 6. ACIS-I light curves (left) and spectra (right) for the seven X-ray sources detected within a 4.4 arcmin radius from the 
center of SI 06, which are bright enough (count rate greater than 2 counts/ks) for a timing and spectral analysis to be carried out. 
The fits to the spectra with absorbed one-temperature plasma model are also shown. 
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Fig. 6. continued - light curves (left) and spectra (right) for the seven brightest X-ray sources in S 106. 



4. Discussion 

4.1. S106 membership 

In the Chandra image of the SI 06 star- forming region 87 X-ray 
sources were detected within a 4.4 arcmin radius of the cen- 
ter of the SI 06 cluster and 6 sources were detected within 0.8 
arcmin radius of the SI 06 south cluster. Of these, 71 sources 
appear to have an IR counterpart in the 2MASS PSC catalogue 
and 22 have no IR counterpart within 3 arcsec from the X-ray 
position. 

The 22 sources without IR counterpart are likely stars, 
members of the SI 06 cluster At a flux level of the order of 
10"''* erg cm"^ s"' the expected number density of background 
sources determined on the basis of the log A^-log S relationship 



for X-ray sources (see e.g. iHasinger et alJ 120011) is 100-200 
sources per square degree. The area of a circle with 4.4 arcmin 
radius is 0.017 square degree, so that the expected number of 
serendipitous extra-galactic sources for a low absorption field 
is between 2 and 3. In practice given the high absorption value 
in this field it is unlikely that any of these 22 sources is a back- 
ground source (of galactic or extra-galactic origin). The fact 
that none of these 22 sources shows up in the soft X-ray band 
is an other indication that most of them are affected by signifi- 
cant absorbing column density and therefore are unlikely to be 
foreground sources. 

From the sample of 71 X-ray sources with IR counterpart 
we selected a subsample of 61 sources which have (photomet- 
rically derived) high extinction levels (A^(H) within a range of 



[0.5 -6.UJ X 10 cm" ) consistent with the sources being mem- 
bers of the clusters. Of all these 61 sources only two (source 30 
and source SI 06 south number 4) are visible in the X-ray soft 
band. This is another, independent, indication of high column 
densities in front of most the sources in the final sample. Of the 
two sources that were detected in the soft band we were able to 
study source 30 in more detail, finding that it has the character- 
istics of a pre-main sequence star (PMS) in the outer edge of 
the association. Thus, the 61 members of the final sample are 
all likely members of the S 106 star forming region. 

4.2. Age of the association 

For the final sample of 61 sources we studied the X-ray lu- 
minosity as a function of mass. As apparent from Fig. |5] and 
Table|3]there is a clear correlation between X-ray luminosities 
and star masses. The presence of a correlation between X-ray 
luminosity and star r nasses has been established fo r other star 
forming regions (e.g. IPreibisch & Zinneckenl2002l for IC 348, 
iFlaccomio et al. 2003 for the ONC) and it is a consequence of 
the fact that members of star-forming regions are young and 
emit at similar, close-to-saturated levels of X-ray luminosity. 
The presence of such a correlation in our sample therefore indi- 
cates that the stars are young and physically associated. A sam- 
ple of field stars at random distances would not present such a 
correlation. 

For stars in the mass bins 0.5-2.0 Mq we derive a median 
valueofLx - 8x10^' ergs"'. Considering the age evolution of 
the X-ray luminosity of solar-mass stars, and using the values 
tabulated by Micela C2002) for the median Lx of a number of 
star forming regions, open clusters and nearby field stars, one 
finds that a median value of Lx = 8x1 0^"^ erg s" ' is very similar 
to the median luminosities in the ONC, the p Oph star-forming 
region and the a Per cluster. This indicates for the S 106 clusters 
an age lower than 3 x 10^ yr In addition, the median X-ray 
luminosity values derived here for the S 106 star- forming region 
are, for all mass bins apar t one (2.0-3.0 Mq-,), wit hin a factor of 
3 of the values derived by Flaccomio et al. ( 2003) for the ONC, 
consistent with the age of SI 06 being simila r to the one of the 
ONC, as well as with the age estimate of iHodapp & Ravnen 
ill99lh ri-2x10^vr1. 

We note that the median value Lx - 10^*^ erg s"' for stars in 
our sample with masses of 0.5-2.0 Mq is not strongly depen- 
dent on the use of a 2 Myr year theoretical isochrone to esti- 
mate extinction coefficients and masses. We have repeated the 
same p rocedure using a zero age main sequence^ by Siess et al.' 
( I2OO0I) and the median luminosity derived for stars within this 
mass range is essentially unchanged. This is because for stars 
with masses of 0.5-2.0 Mq in the 2 Myr hypothesis the differ- 
ence in mass estimate using a zero age main sequence is only 
a factor ~ 2. The value of Lx is affected by the age assump- 
tion via the extinction coefficient and this estimate is also not 
so sensitive to the use of a 2 million year isochrone or a zero 



' Although this is an incorrect assumption for an associations in 
which stars are nearly coeval it allows the dependency of our results 
on the age assumption to be tested. 



age main sequence. In the first case we derive a median value 
ofA/s- - 1.1 and in the second case A^f - 1.5. 

We caution however that in general the use of a zero age 
main sequence rather than the appropriate isochrone can lead to 
significant over-estimate of the cluster total mass. For instance 
in our case the use of a zero age main sequence leads to 22 
stars in our sample having more than 3 Mq while using a 2 Myr 
isochrone one only finds 5 stars. This may explain the sig nifi- 
cant d ifference in cluster mass estimates by Hodapp & Ra vnen 
(Il99li) who derive a total mass of 140 Mq and LDK02 who de- 
rive 400-600 Mq. LDK02 used a li near fi t to a zero age main 
sequence while Hodapp & Ravnerl ( 1199 11) adopted isochrones 
between 0.1-2.0 Myr 

4.3. S106IRS4 

SI 06 IRS 4 is a massive stellar object of spectral type 07- 
BO and luminos ity of [0.4-1] x 10^ Lq (Gehrzetal. 19^ 
iFelli et alJll984l). The line-of-sight extinction is estim ated at 
Av = 20-30 teir oa et alJ Il979': 'Simon et al.' 'l98■3^. From 
2MASS photometric data we estimated for S106 IRS 4 an ex- 
tinction coefficient A^ = 3.35, corresponding to Ay - 30. 
The infrared images are consistent with a model in which 
the recently formed massive stars excites the biconical neb- 
ula from within a large irregular disk of gas and dust which is 
the remnant of the stellar collapse process fGehrz et al. 1 9821 
Harve v et al. 1987) . The near infrared flux from S106 IRS 4 is 
probably not photospheric but results from scattered and ther- 
mal radiation from the inner region of a circumstellar shell. 

Radio observations show the existence of a compact shell 
of radio emission surrounding SI 06 IRS 4 (Sim on et al.lll983t 

II 1 o 

Kurtz et al . 1994). This is interpreted as emission from a flow- 
ing ionized envelope from a short-lifetim e-phase of the pre- 
main sequence evolution of a massive star JSimon et alJll983t 

I — '■ — ni 1 

IFelli et al.MI9841) - an evolutionary state with high mass loss 
which follows an earlier accretion phase. 

The rate at which SI 06 IRS 4 is losing mass is atypi- 
cally high for a star of this bolometric luminosity. At M ~ 
1.6 X 1 0"*" Mfn yr y ' , corresp onding to M/L ^ [1.6-8] x 
1O"""M0 Lq' yr-' dFelli et all 1984). this is 1-2 orders of mag- 
nitude higher than most normal early type stars and comparable 
to the values measured in Wolf-Rayet stars. At the same time 
the terminal wind velocity of Vco ~ 200 km s"' is lower than 
typical values measured f or equally luminous stars, fo r which 
Vco ~ 1000-1500 km s' JPanagia & Macc hettolll982l) . Radio 
obs ervations suggests that thi s wind may be mainly equato- 
rial dHoare et al.ll994 . ISchneider et al.l(l2002l) attribute the dy- 
namics of the molecular gas in the S 106 region to the impact of 
the ionized wind of SI 06 IRS 4, driving a shock into an inho- 
mogeneous molecular cloud. 

The Chandra observation provides the first detection of 
S106 IRS 4 in X-rays. We estimated for this source an X-ray 
luminosity of 2 x lO-'^erg s"' assuming a plasma temperature 
kT = 2.16 keV, as for all the other sources in our sample. This 
luminosity corresponds to Lx/Lboi ~ [0.5 - 1] x 10"*^, which is 
about one order of magnitude lower than typically observed for 
older massive stars JSciortino et alJll99Qi) . although there is a 



large scatter (±1 dex) on the Lx/Lboi relation tor massive stars 
(e.g. Moffat et al. 2002) . The assumed X-ray plasma tempera- 
ture for SI 06 IRS 4 is howeve r likely too high for an early type 
star Chlebo wski et al.l Il989) have shown that typical plasma 
temperatures for wind related emission in O stars are around 
kT - 0.5 keV, significantly lower than the value we have as- 
sumed. Given the high absorption toward the source, the choice 
of plasma temperature is rather critical for determining the X- 
ray luminosity. Assuming kT - 0.5 keV one derives for 
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[1.2- 



SI 06 IRS 4 an intrinsic X-ray luminosity of Lx 
10^' erg s"', corresponding to a range Lx/^boi 
31 X 10~^, at the low er end of the typical values for O stars 
llSciortino e|:^lJ| 1 ^901). On the b asis of the correlations deter- 
mined bv lSciortino et al.l ( 1 19901) in their study of X-ray emis- 
sion from O-stars an X-ray luminosity of 4.8 x 10^' erg s ' is 
a factor of 15 below the value typical for stars with compara- 
ble mass loss, but it is similar to the X-ray luminosities found 
in relation to the star wind momentum flux {F^ = Mvoo = 
2 x 10^^ g cm s~2). 

iKohno et al.l ll2002» have recently reported the detection 
of X-ray emission from four high-mass YSOs in Mon R2, 
deriving typical best fit plasma temperatures, absorption col- 
umn densities and X-ray luminosities of ~ 2 keV, ~ 5- 
10 X 10^^ cm"^ and lO-^'-lO-^'erg s"', that is X-ray luminosi- 
ties similar to the one derived here for SI 06 IRS 4. The X-ray 
flux from the Mon R2 high mass YSOs appears to be highly 
variable with flare-l ike beha vior; because of this and the high 
plasma temperatures lKohno et al. ( 2002) suggest that the X-ray 
activity of these massive YSOs may be magnetically driven, in 
a similar way to that seen in low mass PMS stars. 

Unfortunately, for SI 06 IRS 4 no plasma temperature can 
be determined, which means that we cannot compare it with the 
values derived for the Mon R2 sources and that we have to base 
the estimate of its X-ray luminosity on an assumed plasma tem- 
perature. Nevertheless, assuming a plasma temperature typical 
for older massive stars, the X-ray luminosity of SI 06 IRS 4 is 
consistent to the values predicted for older stars on the basis 
of their wind momentum flux, which is the dominating fac- 
tor in determining the X-ray luminosity of a massive star (as 
established bv .Sciortino et al...l990. - Fig. 16a). This suggests 
that the activity in SI 06 IRS 4 maybe wind-driven, with no 
need to invoke the presence of magnetically confined plasma. 
In SI 06 IRS 4, we would thus be witnessing the onset of X-ray 
emission from the wind, at a stage in which the protostar is still 
deeply embedded in the circumstellar material. 



4.4. Other interesting sources 

Among the 7 sources for which we were able to study the X- 
ray spectra, sources 32, 39 and 60 are intermediate mass stars, 
with estimated masses of 2.4, 2.7 and 2.9 Mq, respectively. 
Were they main sequence they would correspond to A stars, 
which are, at most, weak X-ray sources (Lx ^ 3 x 10^^ erg s"', 
see e.g. pavata & Micela 2003). Given their young age how- 
ever they could be Herbig Ae/Be stars, which have significant 
X-ray activity, or even their precursors, since according to the- 



oretical models a two million year old star ot I-i Mq will have 
spectral type K-G. 

The luminosities (Lx = 2-9 x 10""^ erg s"') and plasma 
temperatures (kT ~ 2 keV) of sources 32, 39 and 60 
are similar to those typically found for Herbig Ae/B e stars 
JPreibisch & Zinnecker 1996: iHamaguchi et alJl2002l) . While 
these values are not inconsistent with the X-ray emission being 
from unseen low-mass counterparts, evidence is gathering that 
the ori gin of the emission are the Herbig Ae/Be stars them- 
selves ( Preibisch & Zinnecken Il996t IHamaguchi et al.l l2002t 
iGiardino et al-2004j). 

Recently iBeuther et alJ (120021) have reported X-ray emis- 
sion from four intermediate mass YSOs in the massive star 
forming region IRAS 19410-1-2336 which have comparable X- 
ray luminosity (Lx = lO-^'-lO-'^ erg s"') and high plasma tem- 
perature (kT ^ 2 keV). 



5. Conclusions 

The Chandra X-ray observation, combined with the 2MASS 
data, has allowed the SI 06 and SI 06 south clusters to be stud- 
ied in more detail confirming them as sites of recent star for- 
mation, with age comparable to that of the ONC. In addition, 
the X-ray observation has allowed the low-mass YSO popula- 
tion of SI 06 to be identified, opening the way for the IR study 
of the individual stars. The X-ray characteristics of this pop- 
ulation appear to be similar to the ones of the (much better 
characterized) ONC. 

We have detected X-ray emission from SI 06 IRS 4, a 
highly embedded massive young stellar object with an excep- 
tional wind for its luminosity. While somewhat dependent on 
the plasma temperature (which cannot be determined from our 
data), the X-ray luminosity of SI 06 IRS 4 appears low in com- 
parison with massive main sequence and more evolved mas- 
sive stars, both on the basis of its bolometric luminosity as 
well as its mass loss rate. Nevertheless, when wind momen- 
tum flux is considered, which includes the influence of both 
the mass loss rate and the terminal velocity, the X-ray lumi- 
nosity of S 106 IRS 4 appears typical of more evolved massive 
stars, suggesting the same process to be at work. The lower lu- 
minosity with respect to stars of similar mass loss rate is due to 
the peculiarly low wind terminal velocity. This would therefore 
suggest that the wind-driven X-ray emission of massive stars is 
already present d uring this early pr e-main sequence phase. As 
discussed e.g. bv'Ballv et ai](ll998l) . S106 IRS 4 appears to be 
extremely young, and Fel li et alJ (11984) argue that the star is in 
a short-lived phase of its early evolution which follows an ear- 
lier accretion phase and it is characterized by high mass loss. 

Further observations of SI 06 IRS 4 and other objects of 
this type (the present observation being thus far unique) will 
be required to establish whether such an early start of the X- 
ray emission in massive stars is a common phenomenon and 
whether the emission is wind-driven already at these early 
stages. Were this the case one could speculate that a wind- 
driven X-ray emission in SI 06 IRS 4 is a consequence of it 
being in a more advanced evolutionary p hase in respec t to the 
massive YSOs detected in the X-ray by iKohno et alJ (120021) 



and tor which they argue that a magnetic conhned plasma is 
present. 
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lable 1. A-ray sources detected in the Chandra data within 4.4 arcmin rrom the center or SlUo. ror those sources tor which 
an IR counterpart in the 2MASS PSC exists the J, H, K magnitude are given together with the separation (r) from the possible 
counterparts. The 's' flag indicates sources which are also detected in the soft band (AE = 0.2 - 1.0 keV). The value of log L^ 
derived from Eq. 3 for the sources with IR counterpart is also given. Seven sources are inconsistent with the estimated age and 
distance of S 106, therefore, their luminosity cannot be computed using Eq. 3. We have marked these sources with 'f. 



Source 


RA(J2000) 


Dec(J2000) 


Count rate 

(ks)-' 


J 


H 


K 


r 
(arcsec) 


logLx 
(log ergs"') 


Soft 


1 


20:27:06.1 


37:23:29 


0.59 ±0.18 


14.5 


13.3 


13.0 


0.68 


29.6 




2 


20:27:07.4 


37:21:11 


0.40 ±0.14 


14.4 


13.2 


12.7 


0.63 


29.7 




3 


20:27:08.5 


37:25:06 


0.35 ±0.12 


14.9 


13.3 


12.5 


1.05 


29.9 




4 


20:27:09.0 


37:22:18 


0.37 ±0.13 


12.9 


12.0 


11.6 


2.29 


29.4 




5 


20:27:09.2 


37:23:31 


0.17 ±0.09 


13.3 


12.3 


11.9 


0.40 


29.2 




6 


20:27:09.2 


37:24:32 


0.55 ±0.17 


16.3 


15.1 


14.5 


0.76 


29.7 




7 


20:27:09.6 


37:21:38 


0.25 ±0.18 


13.6 


12.6 


12.2 


0.97 


29.3 




8 


20:27:09.9 


37:22:28 


0.57 ±0.18 


16.1 


14.4 


13.4 


0.44 


30.1 




9 


20:27:10.7 


37:24:42 


0.16 ±0.07 


15.8 


14.5 


13.9 


0.28 


29.2 




10 


20:27:12.1 


37:23:60 


0.44 ±0.15 


15.0 


13.1 


12.0 


0.28 


30.2 




11 


20:27:12.2 


37:22:53 


0.18 ±0.08 


15.5 


14.1 


13.6 


0.54 


29.3 




12 


20:27:12.5 


37:25:02 


0.48 ±0.16 


16.1 


14.5 


13.9 


0.63 


29.9 




13 


20:27:13.1 


37:22:07 


0.16 ±0.08 


17.2 


15.2 


14.4 


0.45 


29.6 




14 


20:27:13.3 


37:22:13 


0.49 ±0.15 


- 


- 


- 


- 


- 




15 


20:27:14.9 


37:22:43 


0.21 ± 0.08 


16.0 


14.2 


13.6 


0.35 


29.6 




16 


20:27:16.4 


37:22:02 


0.28 ±0.12 


15.5 


14.8 


14.6 


0.47 


f 




17 


20:27:18.0 


37:24:39 


1.29 ±0.19 


16.0 


14.0 


13.2 


0.42 


30.6 




18 


20:27:18.3 


37:22:23 


0.34 ±0.14 


- 


- 


- 


- 


- 




19 


20:27:18.4 


37:21:04 


0.20 ± 0.09 


- 


- 


- 


- 


- 




20 


20:27:19.2 


37:22:36 


0.72 ± 0.21 


15.0 


13.5 


12.9 


0.65 


30.1 




21 


20:27:20.8 


37:23:13 


0.23 ±0.11 


15.1 


13.3 


12.4 


0.37 


29.8 




22 


20:27:21.3 


37:23:43 


3.53 ± 0.46 


17.4 


15.1 


13.8 


0.18 


31.1 




23 


20:27:21.4 


37:24:28 


0.23 ± 0.09 


16.8 


14.9 


13.7 


0.64 


29.8 




24 


20:27:22.0 


37:23:53 


0.18 ±0.08 


16.0 


13.8 


12.5 


0.25 


29.8 




25 


20:27:22.5 


37:19:50 


0.78 ±0.15 


- 


- 


- 


- 


- 




26 


20:27:22.8 


37:23:52 


1.48 ± 0.20 


13.4 


12.1 


11.6 


0.22 


30.3 




27 


20:27:22.8 


37:24:14 


1.05 ±0.17 


14.2 


13.6 


13.3 


0.09 


f 


s 


28 


20:27:23.1 


37:23:38 


0.29 ±0.11 


13.8 


12.2 


11.1 


0.15 


29.9 




29 


20:27:23.3 


37:23:26 


0.83 ±0.16 


13.8 


11.1 


9.0 


0.45 


30.8 




30 


20:27:23.8 


37:22:45 


2.12 ±0.23 


13.3 


11.9 


11.6 


0.28 


30.5 


s 


31 


20:27:23.8 


37:22:09 


0.70 ±0.31 


14.1 


12.8 


11.9 


1.22 


30.2 




32 


20:27:24.0 


37:23:07 


2.31 ±0.24 


11.3 


10.6 


10.2 


0.23 


30.2 




33 


20:27:24.4 


37:23:10 


0.49 ±0.16 


15.2 


13.5 


11.5 


0.38 


30.3 




34 


20:27:24.4 


37:23:40 


0.32 ±0.12 


15.4 


14.0 


13.6 


0.22 


29.5 




35 


20:27:24.6 


37:18:29 


0.67 ±0.19 


13.4 


13.0 


12.9 


0.37 


f 


s 


36 


20:27:24.6 


37:23:25 


0.32 ±0.12 


14.0 


12.5 


12.0 


0.08 


29.7 




37 


20:27:25.1 


37:22:48 


0.32 ±0.13 


- 


- 


- 


- 


- 




38 


20:27:25.2 


37:22:51 


0.65 ±0.19 


14.0 


13.0 


12.9 


0.28 


28.8 




39 


20:27:25.2 


37:23:14 


2.18 ±0.24 


14.3 


12.3 


11.2 


0.10 


30.9 




40 


20:27:25.7 


37:22:57 


0.40 ±0.13 


- 


- 


- 


- 


- 




41 


20:27:26.1 


37:22:59 


0.41 ±0.13 


13.0 


12.0 


11.0 


2.06 


29.9 




42 


20:27:26.1 


37:22:42 


0.17 ±0.08 


- 


- 


- 


- 


- 




43 


20:27:26.2 


37:22:32 


0.56 ±0.17 


13.5 


11.1 


10.3 


0.74 


30.4 




44 


20:27:26.3 


37:22:49 


0.19 ±0.11 


- 


- 


- 


- 


- 




45 


20:27:26.3 


37:22:47 


0.76 ± 0.21 


- 


- 


- 


- 


- 




46 


20:27:26.3 


37:23:31 


0.26 ±0.10 


15.1 


13.2 


12.2 


0.61 


29.9 




47 


20:27:26.5 


37:24:30 


0.41 ±0.13 


15.7 


14.1 


13.6 


0.57 


29.8 




48 


20:27:26.5 


37:22:51 


0.34 ±0.12 


- 


- 


- 


- 


- 




49 


20:27:26.5 


37:23:05 


0.36 ±0.14 


12.3 


10.8 


11.0 


1.19 


29.5 




50 


20:27:26.8 


37:22:48 


0.30 ±0.11 


10.4 


7.7 


5.9 


0.44 


30.3 





lable 1. (continued) X-ray sources within 4.4 arcmin or S1U6 center 



Source 


RA(J2000) 


Dec(J2000) 


Count rate 

(ks)-' 


J 


H 


K 


r 
(arcsec) 


logLx 
(log ergs-') 


Soft 


51 


20:27:26.8 


37:22:43 


0.30 ±0.11 


- 


- 


- 


- 


- 




52 


20:27:27.0 


37:22:53 


0.23 ± 0.09 


- 


- 


- 


- 


- 




53 


20:27:27.0 


37:23:16 


0.27 ±0.11 


- 


- 


- 


- 


- 




54 


20:27:27.1 


37:22:55 


0.98 ±0.16 


12.2 


11.9 


10.5 


2.08 


30.2 




55 


20:27:27.1 


37:23:23 


0.43 ±0.14 


13.9 


12.3 


11.6 


0.49 


30.0 




56 


20:27:27.1 


37:22:45 


0.19 ±0.08 


- 


- 


- 


- 


- 




57 


20:27:27.1 


37:22:36 


0.53 ±0.16 


- 


- 


- 


- 


- 




58 


20:27:27.2 


37:23:23 


0.39 ±0.13 


13.9 


12.3 


11.6 


1.54 


29.9 




59 


20:27:27.2 


37:23:02 


0.26 ±0.11 


13.5 


11.0 


9.8 


2.18 


30.1 




60 


20:27:27.6 


37:22:43 


2.45 ± 0.25 


14.0 


12.6 


10.8 


3.02 


31.0 




61 


20:27:27.7 


37:22:34 


0.35 ±0.13 


- 


- 


- 


- 


- 




62 


20:27:27.9 


37:22:36 


0.42 ±0.14 


14.6 


13.1 


9.5 


2.12 


30.5 




63 


20:27:28.0 


37:21:29 


0.21 ± 0.08 


- 


- 


- 


- 


- 




64 


20:27:28.0 


37:22:53 


0.43 ±0.15 


- 


- 


- 


- 


- 




65 


20:27:28.2 


37:24:47 


0.96 ±0.15 


164 


14.6 


13.7 


0.16 


30.3 




66 


20:27:28.5 


37:24:03 


0.28 ±0.10 


14.8 


13.2 


12.6 


0.35 


29.7 




67 


20:27:28.8 


37:19:50 


0.68 ± 0.20 


134 


12.8 


12.6 


0.21 


f 


s 


68 


20:27:28.8 


37:19:45 


3.71 ±0.30 


13.7 


13.2 


12.8 


0.39 


f 


s 


69 


20:27:28.8 


37:22:42 


0.44 ±0.15 


15.0 


13.4 


12.8 


0.71 


29.9 




70 


20:27:28.9 


37:23:01 


0.62 ±0.11 


- 


- 


- 


- 


- 




71 


20:27:29.2 


37:23:15 


0.27 ±0.10 


- 


- 


- 


- 


- 




72 


20:27:29.2 


37:22:46 


2.48 ± 0.25 


14.7 


12.9 


12.1 


0.10 


30.8 




73 


20:27:29.3 


37:23:19 


0.20 ± 0.08 


154 


13.9 


12.6 


0.50 


29.8 




74 


20:27:29.5 


37:23:40 


0.98 ±0.16 


14.9 


13.1 


12.2 


0.25 


30.4 




75 


20:27:29.8 


37:22:46 


0.28 ±0.10 


15.6 


13.9 


13.0 


0.29 


29.8 




76 


20:27:29.8 


37:22:57 


0.17 ±0.09 


- 


- 


- 


- 


- 




77 


20:27:30.2 


37:22:56 


0.30 ±0.11 


- 


- 


- 


- 


- 




78 


20:27:30.3 


37:22:34 


0.71 ±0.25 


14.8 


13.1 


12.2 


0.31 


30.3 




79 


20:27:30.4 


37:22:39 


0.34 ±0.14 


15.5 


13.9 


13.5 


0.21 


29.6 




80 


20:27:30.7 


37:24:29 


0.31 ±0.11 


12.9 


12.4 


12.2 


0.22 


f 


s 


81 


20:27:30.8 


37:24:11 


0.18 ±0.08 


16.3 


14.3 


13.2 


0.57 


29.8 




82-a 


20:27:30.9 


37:23:21 


0.34 ±0.12 


15.5 


14.0 


14.3 


0.21 


29.6 




82-b 


idem 


idem 


idem 


15.6 


14.0 


13.7 


2.63 


28.5 




83 


20:27:31.6 


37:23:08 


0.67 ±0.13 


14.5 


12.9 


12.2 


0.17 


30.2 




84 


20:27:35.5 


37:20:24 


0.36 ±0.13 


11.7 


11.2 


11.1 


2.58 


f 


s 


85 


20:27:40.5 


37:21:47 


0.27 ±0.11 


16.2 


15.1 


14.6 


0.41 


29.2 




86 


20:27:41.6 


37:25:35 


0.29 ±0.12 


- 


- 


- 


- 


- 




87 


20:27:42.7 


37:23:40 


1.06 ±0.17 


14.0 


13.2 


12.8 


0.46 


29.2 


s 



Table 2. X-ray sources detected in the Chandra data using Pwdetect within 0.8 arcmin from the center of SI 06 south. For those 
sources for which an IR counterpart in the 2MASS PSC exists the J, H, K magnitude are given. The 's' flag indicates sources 
which are also detected in the soft band (AE - 0.2 - 1 .0 keV). The value of log Lx derived from Eq. 3 is also given. 



Source 


RA(J2000) 


Dec(J2000) 


Count rate 

(ks)-' 


J 


H 


K 


r 
(arcsec) 


logLx 
(log ergs-') 


Soft 


1 


20:27:20.0 


37 


17:16 


0.51 ±0.17 


14.6 


13.5 


13.1 


0.57 


29.5 




2 


20:27:22.8 


37 


17:55 


0.64 ±0.19 


15.2 


13.2 


11.6 


0.27 


30.4 




3 


20:27:24.4 


37 


17:41 


0.38 ±0.15 


15.8 


14.1 


13.6 


0.79 


29.8 




4 


20:27:24.8 


37 


17:33 


1.28 ± 0.20 


11.0 


10.2 


9.8 


0.58 


30.1 




5 


20:27:25.4 


37 


17:08 


0.82 ± 0.23 


16.5 


14.3 


12.6 


1.71 


30.6 


s 


6 


20:27:26.2 


37 


17:33 


0.50 ±0.17 


14.8 


12.7 


11.4 


0.44 


30.3 





lable 4. Summary or the derived X-ray properties tor the 7 brightest sources in our sample. A'(H) phot, gives the absorbing 
column density as derived from the IR photometric data (see Sect. l3.2l for the procedure). A^(H)and kT give the best fit values for 
the absorbing column density and the plasma temperature as derived from spectral fitting of the ACIS-I data with one-temperature 
plasma models. The models' reduced ;t'^ and null hypothesis probability are given in the two next columns. F^ and "Fx unabs." 
give the observed and unabsorbed X-ray flux for the best fit model. L]^ and L^ are, respectively, the luminosity derived from 
"Fx unabs" and the luminosity derived with the same approach as for all the other sources in our sample for which a spectral 
analysis was not possible. M gives the stars' estimated mass, and K-S indicates the probability of constancy according to the 
Kolmogorov-Smirnov statistical test. Unless otherwise specified, fluxes are in the energy range 0.8-7.5 keV. A^(H) is in units of 
N22 - 10^^ cm"^, fluxes are in unit of F-o = 10"'-' erg cm"- s"' and Lx in unit of L30 = lO-'^erg s"'. 



Source 


A'(H) phot. 


MH) 


kT 


x' 


P 


F. 


F^ unabs. 


M 


^; 


4 


K-S 




N22 


N22 


keV 






F-13 


F-u 


Mq 


^30 


-^^30 




22 


3.3 


9.7±5.5 


2.64±2.46 


0.62 


0.79 


0.7° 


2.6" 


0.6 


11 


13 


0.25x10-'' 


30 


1.3 


0.8±0.4 


6.71±5.57 


1.00 


0.46 


0.3 


0.3 


1.4 


1.3 


3.1 


0.20x10-2 


32 


0.6 


0.9±0.4 


2.61±1.02 


0.66 


0.87 


0.3 


0.4 


2.4 


1.7 


1.6 


0.22 


39 


3.1 


4.1±1.0 


1.23±0.37 


1.20 


0.25 


0.3 


2.0 


2.7 


8.6 


7.6 


0.56x10-2 


60 


3.3 


3.2±0.8 


2.33±0.75 


0.72 


0.82 


0.4 


1.0 


2.9 


4.3 


9.2 


0.83x10-' 


68* 


- 


0.2±0.1 


1.02±0.11 


1.14 


0.19 


0.2 


0.3 


- 


0.4 


- 


0.59x10-'' 


72 


2.4 


2.1±0.6 


2.99±1.27 


0.79 


0.73 


0.4 


0.8 


1.5 


3.5 


6.7 


0.18 



"range 1.4-7.5 keV 

* likely a foreground source: it is inconsistent with a 2 Myr isochrone at 600 pc 



